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Is the Exotic Hadron X(3872) a D^D Molecule: Precision Determination of the 

Binding Energy of X(3872) 
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It has been proposed that the recently discovered archetypical "exotic" meson, X(3872), with 
A/(X(3872)) = 3871.68 ± 0.17 MeV/c^ and an extremely narrow width, r(X(3872)) < 1.2 MeV, is 
a hadronic molecule of bound D° and D*" mesons. If true, this would establish a new species of 
hadrons, distinct from qq mesons and qqq baryons. It is put to an important experimental test by 
making a high precision measurement of the proposed molecule's binding energy. Using 818 pb~^ of 
e'^e~ annihilation data taken with the CLEO-c detector at 7/;(3770), the decays KsK'^ K~ and 

D^(lf) — > K^TT^-K^TT~ have been studied to make the highest precision measurement of D'^ mass, 
M{D") = 1864.851±0.020±0.019±0.054 MeV/c^ where the first error is statistical, the second error 
is systematic, and the third error is due to uncertainty in kaon masses, or 1864.851 ±0.061 MeV/c^ 
with all errors added in quadrature. This leads to M(D° + Z)*°) = 3871.822 ± 0.140 MeV/c^, 
and the binding energy BE(X(3872)) = M{D° + D'°) ~ M(X(3872)) = +142 ± 220 keV. At the 
90% confidence level this leads to the conclusion that X(3872) is either unbound by as much as 
140 keV, or it is bound by less than 420 keV. If bound, X(3872) has a very large radius; the central 
value of binding energy corresponds to a radius of 12 fm, and the lower limit to 7 fm, both being 
uncomfortably large for a molecule. 

PACS numbers: 14.40.Lb, 12.40.Yx, 13. 25. Ft 

Recent observations at the B-factories of many un- 
expected resonances, loosely called "exotic", have given 
rise to great excitement in the spectroscopy of heavy- 
quark hadrons. In the mass region above bound char- 
monium resonances, M > 3.73 GeV/c^, most of the ob- 
served resonances have widths which range from 30 to 
200 MeV, and many remain potential charmonium candi- 
dates. However, one resonance, dubbed X(3872), has ac- 
quired the status of the quinessential exotic because of its 
unique properties. It has been observed in many diverse 
experiments, by Belle P, i, BaBar 0,13, CDF [Hi], 
D0 [7], LHCb d, and CMS % The closeness of its mass, 
Af(X(3872)) = 3871.68 ± 0.17 MeV/c^ (T^ to the sum 
of the masses of two open charm mesons D^{J^ ~ 1^) 
and D*^{J^ = 1^), and its extremely narrow width, 
r(X(3872)) < 1.2 MeV (90% CL), has given rise to the 
proposal that it is a D^D*^ molecule. Many of its decays 
have been measured, and it is found that all decay final 
states invariably contain a charm quark and an anticharm 
quark, which would suggest that it is a narrow charmo- 
nium resonance. Angular correlation measurements limit 
its J^<^ to 1++ or 2-+ [nl m, so that the likely char- 
monium states would be 2'^Pi(l++)x^]^ or \^D2{2 ^)rjc2- 
Unfortunately, the predicted masses of both these states 
are quite far (- -h75 MeV/c^, and 40 MeV/c^, re- 
spectively) from 3872 MeV/c^, which makes it difficult 
to identify X(3872) as a pure charmonium state. 

Numerous theoretical models for X(3872) have been 
proposed, and several reviews of the different possibili- 
ties exist in the literature However, despite some 
problems, the most popular explanation for X(3872) re- 



mains that it is a loosely bound molecule of the and 
D*^ mesons. If this explanation is correct, X(3872) would 
be member of a new species of hadrons, distinct from qq 
mesons and qqq baryons. This indeed would be a most 
dramatic development in hadron spectroscopy, one that 
needs to be submitted to critical scrutiny. 

Obviously, one of the most important properties of a 
molecule is its binding energy and it is necessary to make 
an accurate determination of it. Since the difference be- 
tween the masses of the and D*° mesons has been 
accurately measured to be 142.12 ± 0.07 MeV/c^ [l3|, a 
precision determination of the binding energy of X(3872) 

as a D molecule requires the highest precision mea- 
surement of the mass of the 13° meson. In this letter we 
report on such a measurement. 

We had earlier 14] reported the 
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measurement 
Ks(j), (f> K+K- us- 
ing 280 pb"^ of CLEO-c data taken at the ■0(3770). 
We determined M{D°) = 1864.847 ± 0.178 MeV/c^, 
which led to an uncertainty of ±363 keV/c^ in the mass 
[M(D°) + M(D*°)]. With the then known value of 
M(X(3872)) which had an uncertainty of ±500 keV/c^, 
the binding energy was determined to be BE(X(3872)) = 
M{D° + D*°) - M(X(3872)) = 600 ± 600 keV/c^. Since 
then, several improved measurements of M(X(3872)) 
have been made [2, 0, S, H @| , with the present PDG av- 
erage M(X(3872)) = 3871.68 ± 0.17 MeV/c^ [T^, and it 
is now necessary to make a correspondingly more precise 
measurement of [M{D'^) + M{D )] in order to determine 
BE(X(3872)) with higher precision. In this letter, we re- 
port such a measurement, which raises serious questions 
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for the \D°D*°) molecule model of X(3872). 

A nearly factor three improvement in the precision of 
M{D^) has become possible because of two reasons. We 
now have nearly three times more CLEO-c data avail- 
able, 818 pb~i of data taken at the ■0(3770), y/s = 
3770 MeV, and in addition to D° KsK+R-, we study 
the nearly forty times more prolific decay, D'^{D ) — >■ 
K^TT^TT+TT- . The data taken at = 3770 MeV is ide- 
ally suited for this measurement because "0(3770) decays 
almost exclusively to DD (branching fraction — 93lg%), 
and the D-mesons are produced almost at rest. 

In the present investigation we determine the mass of 
with a precision of ~ 60 keV/c^. This requires im- 
provement of the default solenoid magnetic field calibra- 
tion of CLEO-c, and to track its small variation with 
time. We do the calibration by choosing to anchor our 
mass measurements to the high precision measurements 
of the masses of 0(25') and J /ip with uncertainties of 
±15 keV/c^ and ±12 keV/c^, respectively, made by the 
KEDR Collaboration by the resonance-depolarization 
technique [isl [l6| . Our investigation involves several 
steps. We first recalibrate the CLEO-c solenoid mag- 
netic field using the KEDR masses in a study of the ex- 
clusive decay ■0(25') tt+tt" J/^0 using CLEO-c data for 
25 million ■0(25'). With the recalibrated magnetic field 
we make a precision measurement of the mass of Ks in 
the inclusive decay, 0(25) Ks + X. Using M{Ks) 
so determined we do fine tuning of the magnetic field 
for each individual CLEO-c dataset at -^s = 3770 MeV 
via the inclusive decay D — > Ks ± X. We use these 
fine tuned fields to make our measurements of mass 
in the two exclusive decays: — > KsK^K^ and 

The data were taken with CLEO-c detector [17;] , which 
consists of a CsI(Tl) electromagnetic calorimeter, an in- 
ner vertex drift chamber, a central drift chamber, and 
a ring imaging Cherenkov (RICH) detector, all inside a 
superconducting solenoid magnet providing a 1.0 Tesla 
magnetic field. For the present measurements, the im- 
portant components are the drift chambers, which pro- 
vide a coverage of 93% of 47r for the charged particles. 
The detector response was studied using a GEANT-based 
Monte Carlo (MC) simulation including radiation correc- 
tions [18]. 

For the analysis of 0(25) decays, ip{2S) — ^ tt+tt" J/0, 
J/ip iJ.'^iJ.', and ^/'(25) Ks + X, Ks ^ tt+tt^, 
we select events with well-measured tracks by requir- 
ing that they be fully contained in the barrel region 
(I cos 1 < 0.8) of the detector, and have transverse 
momenta > 120 MeV/c. For the pions from Ks de- 
cay, we make the additional requirement that they orig- 
inate from a common vertex displaced from the inter- 
action point by more than 10 mm. We require a Ks 
fiight distance significance of more than 3 standard de- 
viations. We accept Ks candidates with mass in the 



range 497.7 ± 12.0 MeV/c^. We identify muons from 
J/ip decays as having momenta more than 1 GeV, and 
Ecc/p < 0.25 for at least one muon candidate, and 
Ecc/p <0.5 for the other muon. We require that there 
should be only two identified pions and two identified 
muons with opposite charges in the event. The momenta 
of pairs is kinematically fitted to the KEDR J/i/j 

mass, M{J/tp) = 3096.917 MeV/c^, and only events with 
< 20 are accepted. We also require that there should 
not be any isolated shower with energy more than 50 
MeV in the event. 

As stated earlier, to make a precision recalibration of 
the solenoid magnetic field we reconstruct ■0(25) in the 
decay V(25) tt+tt" J/V-, J/V' ^ M+A*". Using KEDR 
masses, M(0;(25))kedr = 3686.114± 0.015 MeV/cf [Til 
and M(J/V')kedr = 3096.917 ± 0.012 MeV/c^ [3, we 
determine the magnetic field correction required to mod- 
ify the pion momenta such that M(^0 (25)) present becomes 
identical to M(^0(25))kedr- With pions with momenta 
< 600 MeV/c, the required correction is determined to be 
+0.029% (or 0.29 Gauss) in the default CLEO calibration 
of the magnetic field. Fig. l(top) shows the distribution 
of AM(V'(25)) = M(V'(25)) 

Present 

)-Af(V'(25)) KEDR af- 
ter the correction. We make fits to the unbinned data 
in the full range with peaks parameterized as sum of a 
simple Gaussian function and a bifurcated Gaussian func- 
tion, and a linear background. The fit has the number 
of events, iV(V'(25)) = 125300 ± 356, and xVdof=1.00, 
and gives AM(V(25)) = 0.0 ± 6.7 keV/c^. An identi- 
cal procedure is used to fit all other mass distributions 
presented in this letter. 

Having corrected the magnetic field, we use it to an- 
alyze the same ip{2S) data set for the inclusive de- 
cay, ^0(25) — ^ Ks + X, Ks — > TT+vr^ for pions in the 
same momentum region, < 600 MeV/c^. The fit to the 
TT+TT^ invariant mass distribution, shown in Fig. 1 (bot- 
tom), has x^/dof=1.06. It leads to the number of Ks, 
N{Ks) = 256859 ± 739, and 

M{Ks) Present 

= 497.600 ± 0.007(stat) MeV/c^ (1) 

We next select candidates in the data taken at 
y/s — 3770 MeV using the standard CLEO D-tagging cri- 
teria, which impose a very loose requirement on the beam 
energy constrained mass, as described in Ref. 
Because the data at y/s = 3770 MeV were taken in sev- 
eral smaller sub-runs, the solenoid magnetic field needs 
to be corrected for possible small variations from sub-run 
to sub-run. We do so by analyzing each sub-run for the 
inclusive decay D — Kg + X, Ks — tt+tt" and requir- 
ing that the field be corrected to bring each M{Ks) to 
the value Af(ifs)cLEO in Eq- 1- These corrections were 
found to be < ±0.030%, consistent with what was found 
for ip{2S). With these corrections in place, individual 
data sets were analyzed for the decays — > KsK~^ K' , 
Ks — TT+TT", and £)°(D°) K^tt^tt^tt^ . In order 
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FIG. 1. Top plot: AM(i/'(25)) = M{tv+tv- J/iP) - 
M{%I}{2S))kedr for the exclusive decays tp{2S) n^-K~ J/ip, 
Jf-ip — >■ iJi^ ijT using the corrected magnetic field. Bottom 
plot: M{-K^-R~) for the inclusive reaction ip{2S) Ks + X, 
Ks — > Tv^-K^ , using the corrected magnetic field. The curves 
show peak fits with a sum of a simple Gaussian function and 
a bifurcated Gaussian function and a linear background. 



to use the same magnetic field calibration for all decay 
channels, we use only those events in which final state 
pions and kaons have momenta < 600 McV/c. Wc se- 
lect well-measured tracks which have specific ionization 
energy loss, dE/dx, in the drift chamber consistent with 
pion or kaon hypothesis within 3 standard deviations. 
For the Ks candidates from the exclusive 1)° decays, we 
perform a mass-constrained (IC) kinematic fit, and ac- 
cept in our final sample Ks with < 20. 

The final mass distributions for the different sub-runs 
were added together, and total data at y/s = 3770 MeV 
for the decays, KsK^K~ (which includes events 

in which K'^K~ form a </) meson) and D^{D ) ^■ 
ii'^7r^7r"'"7r~ {KZir henceforth) are shown in Fig. 2. The 
distributions were fitted with peaks as described before, 
and linear backgrounds. The fits have full widths at half 
maximum of 6.4 MeV for KsK+K', and 8.1 MeV for 



FIG. 2. Invariant mass spectra for the decays (top) -D" — >■ 
KsK+R-, (bottom) D° K'^-n^n+Tv' . Fitted masses with 
statistical errors only are given. 

K3tt, and x^/dof of 1.02 and 1.04, respectively. The re- 
sults of the fits are listed in Table I. Both results for 
M{D^) are seen to be consistent within statistical errors, 
and their average is 

M{D^)present = 1864.851 ± 0.020(stat) MeV/c^. (2) 



TABLE I. Results of fits to the mass distributions. The sys- 
tematic errors arc as listcxl in Tabic II and described in the 
text. 



Decay 


A'' (events) 


M(D"), MeV/c^ 


D" KsK+R- 


1,655 ±43 


1864.871 ± 0.074 ± 0.063 


D° KSiT 


76, 988 ± 388 


1864.848 ± 0.021 ± 0.061 


Average 




1864.851 ± 0.020 ± 0.057 



The systematic errors in M{Ks) and M{D^) were ob- 
tained as follows. 



4 



For M{Ks) measurement, we have corrected the 
magnetic field using KEDR measured M{%p{2S)) and 
Af(J/'0), which have the total errors of ±15 keV/c^ 
and ±12 keV/c^, respectively 15, The change in 

M{Ks) due to the change in magnetic field is factor 
1.46 smaller than the change in M('0(2S')). We there- 
fore assign (±15/1.46)~ ±10.3 keV/c^, and (±12/1.46)- 
±8.2 keV/c^, as the uncertainties in M{Ks) due to the 
uncertainties in M{ip{2S)) and M{J/ip). The variation of 
the fit range by ±2 MeV/c^ yielded ±4 keV/c^ systematic 
error. Changing the fits to the background from polyno- 
mials of order one to polynomials of order two changes 
M{Ks) by < 1 keV/c^. The effect of the possible for- 
mation of ipi"^^) at an energy different from the peak 
was investigated in detail. This shift was estimated to be 
±7 keV/c^, and it contributes ±5 keV/c^ to the system- 
atic error in Ks mass. The sum in quadrature of all the 
above contributions is a total systematic uncertainty of 
±15 keV/c^ in M{Ks). 

Our final result for M{Ks) is thus 

M{Ks) Present = 497.600 ± 0.007 ± 0.015 McV/c^ (3) 

Here, and elsewhere when mentioned separately, the first 
error is statistical and the second error is systematic. 
With statistical and systematic errors added in quadra- 
ture our resuh M(i^s)prosont = 497.600 ± 0.017 MeV/c^, 
is the world's most precise single measurement of M{Ks), 
as illustrated in Fig. 3(top). The PDG 2012 average 
of all previous measurements is M{Ks) — 497.614 ± 
0.022 MeV/c^. 



TABLE II. Systematic errors in M{D°) for ranee of variation 
of different parameters. The two values of the total systematic 
errors correspond to the uncorrelated systematic errors, and 
the correlated systematic errors due to uncertainties in Ks 
and K"^ mass measurements. 



M(Kg), MeV 



Source 


AM(D»)sy.t 


AA/(D°),y,t 


Systematic Error (keV/c^) 


D° KsK+K- 


D° ^ K3tt 


coseUax: 0.8, 0.75 


23 


12 


Pniin (trans): 120, 135 MeV/c 


26 


8 


Pmax (total): 650, 570 MeV/c 


4 


7 


Fit Range ±10 MeV 


3 


15 


Bkgd. Polynom. (1,2 order) 


1 


1 


MC Input/Output 


5 


5 


Total (uncorrelated) 


35 


22 


Error in Ks mass 


47 


55 


Error in K^ mass 


22 


13 


Total (correlated) 


52 


57 



The systematic errors in M{D'^), as determined by 
varying event selection and peak fitting parameters, are 
summarized in Table II. Part of each AM(D'')syst listed 
in the table can be due to changes in statistics when the 
parameter values are changed, but to be conservative we 
assign the full variations in M(Z?") to the systematic er- 
rors. The correlated systematic errors listed in Table II 
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FIG. 3. Present results and previous mass measurements [13], 
M{Ks) (top) andM(D°) (bottom). The horizontal solid lines 
show PDG(2012) average. The dashed lines show PDG(2012) 
error bands llOll. 



arise from uncertainties in the masses of Ks and K^. 
The ±17 MeV/c^ uncertainty in the mass of Ks used to 
fine-tune the solenoid magnetic field leads to the largest 
uncertainties in M{D°), ±47 keV/c^ in D° KsK+K- 
and ±55 keV/c^ in D° KSn. The PDG(2012) mass 
of K^ has an error of ±13 keV/c^ [l^. It leads to ±22 
keV/c^ and ±13 keV/c^ uncertainties in M{D°) for the 
decay modes of D° KsK+K' and D° Kin. 

In the two decay modes the sums in quadrature of the 
uncorrelated systematic uncertainties are 35 keV/c^ and 
22 keV/c^, and the sum of correlated uncertainties due to 
Ks and K^ masses are 52 keV/c^ and 57 keV/c^. Taking 
proper account these uncertainties, the average Z?" mass 
for the two decay modes is: 

M{D°)present = 1864.851 ± 0.020 ± 0.019 MeV/c^. (4) 

There is an additional uncertainty of ±0.054 MeV/c^ due 
to uncertainty in kaon masses. With all uncertainties 
added in quadrature, our present result 



M{D'')pre 



1864.851 ± 0.061 MeV/c^ (5) 



is the world's most precise single measurement of the 
mass of the meson, as illustrated in Fig. 3(bottom). 
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It supercedes our previous result in Ref. lj| which was 
based on part of the data used in the present inves- 
tigation. Our resuh for M{D'^), and the PDG value 
A[M{D*°) - M{D°)] = 142.12 ± 0.07 MeV/c^ lead 
to 

Af(D° + Z3*°) = 3871.822 ±0.140 MeV/c^ and (6) 



for this privilege. This research was supported by the 
U.S. Department of Energy. 



BE(X(3872)) = +142 ± 220 keV, 



(7) 



using Af(X(3872)) = 3871.68 ± 0.17 MeV/c^ 0. At 
90% confidence level this result corresponds to X(3872) 
being unbound by 140 keV, or being bound by at most 
424 keV. 

As is well known, a universal property of a weakly 
bound system of two constituents with reduced mass fi, 
and binding energy BE is that the root-mean-square sep- 
aration of the constituents, or the "radius" of the com- 
posite, is given hy d = l/i/S/^BE. The central value and 
the 90% CL upper limit of binding energy lead to 

BE(X(3872)) = 142 keV, d(X(3872)) = 12 fm 
BEmax(X(3872)) = 424 keV, dmi„(X(3872)) = 7 f m 

(8) 

With the early binding energy estimates of the order of 
1 MeV [li], which corresponds to d(X(3872)) = 4.5 fm, 
the long-range interaction responsible for the binding of 
X(3872) was suggested to be pion exchange. The present 
determination of BE — 142 ± 220 keV corresponds to a 
radius as large as 12 fm, or at least 7 fm (90% C.L.), 
twice as large as the deuteron, and it is difficult to see 
how pion exchange could explain the binding of D*^ and 
into a molecule of this size. We recall that sev- 
eral other observations have also raised questions for the 
molecular model. These include too large a cross sec- 
tion for X(3872) formation at the Tevatron, too large a 
ratio ct(X(3872) ^ 7V'(25))/ct(X(3872) ->■ 7J/V'), and 
the possibility that JPC(X(3872)) = 2"+ flS^. Together 
with our measurement of the uncomfortably large size, 
the explanation of X(3872) as a D^D molecule appears 
to have serious problems. 

This investigation was done using CLEO data, and as 
members of the former CLEO Collaboration we thank it 
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